Defining the effective temperature of a quantum driven system from current-current 

correlation functions 
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We calculate current-current correlation functions and find an expression for the zero-frequency 
noise of multiterminal systems driven by harmonically time-dependent voltages within the Keldysh 
non-equilibrium Green's functions formalism. We also propose a fluctuation-dissipation relation for 
current-current correlation functions to define an effective temperature. We discuss the behavior of 
this temperature and compare it with the local temperature determined by a thermometer and with 
the effective temperature defined from a single-particle fluctuation-dissipation relation. We show 
that for low frequencies all the definitions of the temperature coincide. 
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I. INTRODUCTION 

Over the last years we have witnessed a technological 
trend towards miniaturization of electronic circuits. This 
tendency has been accompanied by a growing research ac- 
tivity focused on achieving a better understanding of the 
mechanisms for heat dissipation and energy flow in meso- 
scopic systems. However, the motivation for the research 
in this area is not only technological, because the very 
fundamental concepts of standard statistical mechanics 
and thermodynamics are put into test when studying 
these systems, even more when the process under consid- 
eration corresponds to an out-of-equilibrium situation. 

Several efforts have been made towards the extension 
of standard thermodynamical concepts to the out-of- 
equilibrium evolution of different systems. Some well- 
known examples are the aging regime of glassy systems, 
sheared glasses, granular materials and colloids*^— A 
very successful achievement in the characterization of 
such nonequilibrium states has been the identification of 
an effective temperature, i.e., a parameter with the same 
properties of the temperature of a system at equilibrium 
that is useful to describe the evolution of nonequilibrium 
systems. For instance, for glassy systems the definition 
of effective temperature was introduced^ by means of a 
generalization of the equilibrium fluctuation-dissipation 
relations (FDR) and the physical meaning of this concept 
was supported by showing that such a temperature would 
coincide with the one measured by a thermometer^ The 
definition of an effective temperature from a FDR was in- 
troduced for quantum glassy systems in Ref. and later 
analyzed in electronic systems-. More recently, these 
temperatures were studied in an Ising chain after a sud- 
den quench £ 

The physics of the mesoscopic scale is ruled by the 
quantum coherence of the particle propagation. This 
originates non-trivial interference mechanisms and sur- 
prising effects. Well known examples are the violation of 
the Fourier's Law in low-dimensional phononic systems^ 
as well as the 2k p oscillations of the local voltage and 
the negative electrical resistance^ 



In the past few years there have been many exper- 
imental attempts to locally characterize the heat flow 
in non-equilibrium systems. For example, Pothier and 
coworkers^ measured the local energy distribution func- 
tion in metallic diffusive wires in a stationary out-of- 
equilibrium situation. More recently, Altimiras and 
coworkers 10 measured the electron energy distribution in 
an integer quantum Hall regime with one of the edge 
channels driven out-of-equilibrium. Chiral heat transport 
has been investigated in the quantum Hall regime using 
micron-scale thermometers^ and later explained with 
the introduction of a local temperature along the edgeJ^ 
The idea of defining a non-equilibrium local temperature 
has been useful to study out-of-equilibrium transport in 
other mesoscopic systems. For example, thermoelectric 
transport has been studied with the aid of the local tem- 
perature determined by an ideal thermometer ^ Also the 
concept of effective temperature has been useful to study 
heat exchange between a nanojunction and its environ- 
ment, which can act as a freezing agent ^ and to study 
mesoscopic superconductors^ Another example is the 
prediction that a superconducting wire can remain in 
superconducting state even in contact with a bath that 
greatly exceeds the critical temperature if the effective lo- 
cal temperature is maintained below the critical valued 
A local temperature has also been defined to character- 
ize the heat transport in molecular devices^ Additional 
studies have been reviewed in Ref. [TH 

In a previous workiS we defined local and effective tem- 
peratures in electronic quantum systems driven out of 
equilibrium by external ac potentials. Examples of such 
systems are quantum dots with ac voltages acting at their 
walls (quantum pumps) 20 and quantum capacitors^ In 
that work we presented two concepts, which are the lo- 
cal and the effective temperatures. The local temper- 
ature was introduced following a procedure inspired in 
a work by Engquist and Anderson^ The idea is to in- 
clude a thermometer in the microscopic description of 
the system. On the other hand the effective tempera- 
ture is defined from a local FDR involving single-particle 
Green's functions. We showed that for low driving fre- 
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quencies both ways of defining the temperature coincide. 
In a more recent work-2 3 - we slightly generalized the def- 
inition of the thermometer to consider the possibility of 
simultaneously sensing the local temperature and the lo- 
cal chemical potential of the sample. We showed that 
the new local temperature determined by this new def- 
inition coincides with the previous one. Even more, we 
showed that such a parameter verified the thermodynam- 
ical properties of a temperature, meaning that its gradi- 
ent signals the direction for heat flow at the contacts. 



The aim of this work is to analyze the role of effective 
temperatures within the context of a FDR for current- 
current correlation functions. The motivation is twofold. 
On one hand we are interested in testing the robustness 
of the definition of an effective temperature from a FDR, 
at the level of a correlation function different from the 
one we have considered in our previous work. On the 
other hand, current-current correlation functions are par- 
ticularly appealing quantities since they are related to 
noise, which can be experimentally measured and con- 
tain valuable information on the nature of the elementary 
particles that take part in the transport process. The 
zero-frequency noise is usually used to characterize the 
correlations between particles in mesoscopic systems.™ 
Additionally in quantum pumps, noise is related to the 
possibility of having quantized pumping^ and it con- 
tains information that cannot be extracted from the time- 
averaged current.- 6 Current correlations in mesoscopic 
coherent conductors were first discussed by M. Biittiker 
in Ref. [27] and since then an extensive theoretical litera- 
ture on noise in mesoscopic systems analyzed within the 
scattering matrix formalism has been developed] 26 ' 28 " — 
We use here another approach, which is based Keldysh 
formalism. For non-interacting systems both treatments 
were proved to coincide at the level of the description 
of the current for do 3 -! and ac-driven systems.— In the 
present work we show that this is also the case for the 
current fluctuations correlations. The main goal of this 
work is to show that the effective temperature obtained 
from a fluctuation-dissipation relation for current-current 
correlation functions coincides with the local tempera- 
ture defined using a thermometer and thus verifies the 
same thermodynamical properties of the latter. 



This paper is organized as follows. In Sec. HJ we 
present the model and summarize the theoretical treat- 
ment. In Sec. Mil we review three definitions of tempera- 
ture addressed in recent work o 19 ' 23 . In Sec. II VI we derive 
general expressions for current-current correlation func- 
tions and an explicit expression for the zero-frequency 
noise within the Keldysh Green's functions formalism. In 
Sec. |V]we present numerical results for a particular sys- 
tem. Section IVT1 is devoted to discussion and conclusions. 
We give some details of the calculation in the Appendix. 



II. MODEL AND THEORETICAL TREATMENT 

In Fig. Q] we display the same setup as in Refs. [l^ 
and representing a quantum driven system, with the 
Hamiltonian H sys (i) , connected to a probe characterized 
by Hp. The total system is then described by 



H(t) = H sys (t) + H c p + Hp, 



(1) 



with H c p implementing the local coupling between the 
system and the probe. The Hamiltonian corresponding 
to the driven system can in turn be written as 



H sys (t) =H L + H cL + H c (t) + H cI 



Hr 



(2) 



where Hc(t), Hl and Hr stand for the Hamiltonians 
of the central part and left and right reservoirs, coupled 
among themselves via the Hamiltonians H c l and H c p. 

The Hamiltonian describing the central system (C) 
contains the ac fields and can be written as Hc(t) = 
Ho+Hv(t). We assume that Hq is a Hamiltonian for non- 
interacting electrons while Hy(t) is harmonically time 
dependent with a fundamental driving frequency Qq- We 
leave further details of the model undetermined as much 
of the coming discussion is model independent. 

All three reservoirs (left, right and the probe) are mod- 
eled by systems of non-interacting electrons with many 
degrees of freedom, i.e., H a = J2k a £ kacl a c ka , where 
a = L, R, P. The corresponding contacts are described 
by H ca = w ca J2 ka ( c ka c ia + cj a c ka ), where la denotes 
the coordinate of C where the reservoir a is connected. 
As in previous works ; 8 ' 19 ' 23 ' 33 ' we consider non-invasive 
probe and we treat w c p at the lowest order of perturba- 
tion theory when necessary. 

We will analyze the out-of-equilibrium dynamics of 
this system within the Schwinger-Keldysh Green's func- 
tions formalism. Within this formalism, instead of the 
usual time-ordering operator used in equilibrium theory 
a contour-ordering operator which orders time-labels ac- 
cording to their order on the Keldysh contour is intro- 
duced. The single particle propagator reads 



iG jJ ,{t,t l ) = (Tc[c j {t)cUt')\). 



(3) 



The contour-ordered Green's function contains four dif- 
ferent functions depending on where the times t and t' are 
over the Keldysh contour^ It is easy to see that they are 
not all independent. We then consider the lesser, greater 
and retarded Green's functions, 



iGf tj ,(t,t>) 



-<4(*>i(*)>> 

<9(*)ct,(i')), 

Q{t-t'){\c 3 {t),cUt') 



(4) 



where 



denote the anticommutator of the fermionic 



operators, (...) is the quantum statistical average and the 
indexes j, j' denote spatial coordinates of the system. 
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FIG. 1: Scheme of the setup. The central device is a wire 
with two barriers of height Eb connected by its ends to two 
reservoirs (L and R). The third reservoir (P) represents the 
probe, which consists of a macroscopic system weakly coupled 
to a given point of the central device. In this setup, transport 
is induced by two oscillating ac fields (both with the same 
amplitude and frequency but with a phase lag) applied at 
the points where the barriers are located. The Left and Right 
regions depicted in this scheme are related to the heat current 
that flows into the respective reservoirs . 23 ' 33 



These Green's functions can be evaluated after solving 
the Dyson equations. 

In this work we will focus on current-current correla- 
tion functions. The current in reservoir a at time t is 
defined by the operator 32 

J a (t) = iw ca E (4a (*)^a (*) ~ 4a (*)Cfca (tj) , (5) 
ka 

which obeys bosonic commutation rules. The ensuing 
connected contour-ordered propagator reads in this case 

iC a p(t,t') = (Te[J a (t)J ff (t')]) - (J a (t))(Jp(t')), (6) 

while the lesser, greater and retarded Green's functions 
are 

iC> (t,t') = {J a (t)J fj (t'))-(J a (t))(Mt')), 
iC^(t,t') = Q(t-t'){[j a (t),Mt')\_). (7) 

where [,]_ denote the commutator of the currents. 

For the case of harmonic driving it is convenient to 
use the Floquet-Fourier representation of the Green's 
functions^ 



Ajji (t, t 



E / ^e-^+^A oA k^). 

J — OO 



k=—oo 



where A stands for single-particle © or current-current 
© propagators. 



In general the Keldysh and retarded Green's func- 
tions, can be expressed in terms of the lesser and greater 
Green's functions via 

A*,(M') = A> r (t,t') + A< 3 ,(t,t'), 

A%, (t, = 9(t - [Af j, (t, - A< ., (t, t')] . (9) 

From the definition given in Eq. ([5]) it is straightfor- 
ward to see that the Floquet-Fourier components of these 
functions can be written as 



duj> Al f (k,uj') - A< f (k,uj') 



2tt 



uj - uj' + i0+ 



(10) 



III. DEFINING THE TEMPERATURE 



A. Local temperature determined by a probe 

In Ref. [l^ we defined the local temperature (Tip) of 
the site IP of the system as the value of the temperature 
of the probe such that the time-averaged heat exchange 
between the central system and the probe vanishes. 

It can be shown^ that, given Hc(t) without many- 
body interactions, the dc component of the heat current 
flowing from the central system to the thermometer can 
be expressed as(fi=fcs = e = l) 



J 



Q 



E E 



dtu 



27r {[/» 



a.—L,R,P k— — oo ' 

x( Wfc - fi)r P (uj k )T a (uj) \G? PM (k, }, (11) 

where u>k = uj + kilg, while T a (oj) — — 2ir\w a \ 2 ^2 ka 5(to — 
Ska) are the spectral functions characterizing the reser- 
voirs {a = L,R,P), and f a (ui) = l/feM"-^) + X ] is 
the Fermi function, which depends on T„ = l//3 a and 
Ha respectively the temperature and the chemical poten- 
tial of the reservoir a. Thus, the local temperature Tjp 
corresponds to the solution of the equation 



J^(T lP ) = 0. 



(12) 



In general, Eq. (|12p must be solved numerically, but 
under certain conditions, an analytical expression can 
be found. In particular, for the low temperature weak- 
driving adiabatic regime, which corresponds to small am- 
plitudes and frequencies of the driving potential, and for 
n < T^. 



'ip 



T 1 + Ag? 



(/x)O 



(13) 



where 

k (n) 



i 



(fc) 



n+2 



d n {w(k,u)} 
duj n 



(14) 
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E 

y=L,R 



Gfi a (fc,w)| T a (oj). 



(15) IV. CURRENT-CURRENT CORRELATION 

FUNCTIONS AND EFFECTIVE TEMPERATURE 



An alternative definition of local temperature was dis- 
cussed in Ref. |23|, where the fact that the heat current 
is related to the charge current was taken into account. 
Then, the local temperature (T* p ) and the local chemical 
potential (fJ-f P ) were defined from the condition of simul- 
taneously vanishing of the time-averaged charge and heat 
currents between the probe and the system. That is 



Jp(Ti*p, Hip) 



0. 
0. 



(16) 



where (see Refs. HSH) 



due 



J p= E E / 

a=L,R,P k=-oo J ~°° 



xr P K)r Q ( w )|Gfp !iQ (fc, w )| 2 }, 



(17) 



is the dc component of the charge current flowing through 
the contact between the system and the probe. 

The simultaneous equations given in Eq. (|16[) can be 
solved numerically for any situation, but an analytical ex- 
pression can be found within the low temperature weak- 
driving adiabatic regime, when Oq T, and leads to 
T* p = T lP , given in Eq. CSJ). 



B. Effective temperature from a single-particle 
fluctuation-dissipation relation (FDR) 



For systems in equilibrium, the fluctuation-dissipation 
theorem establishes a relation between the Keldysh (cor- 
relation) and the retarded Green's functions. In Ref. [l9|, 
we defined a local FDR involving single-particle Green's 
functions from which an effective temperature for the site 
I (Ti ff = l/(3f ff ) can be extracted, 



iGfi{Q,u)) -iGf { (0,^) = tanh 



(18) 



with ipi(uj) = -2Im[G/ i ,(0, w)} = Sfc<M fc , w -fc)- In gen- 
eral, Eq. (TT5|) defines an effective temperature that might 
depend on uj, so the limit u) — > fi is taken. An extra term 
is added to the lhs of Eq. (fT8]) because the rhs is always 
zero at 10 = /i but G^(0, jii) is not necessarily zero in an 
arbitrary out-of-equilibrium situation. 

Within the low temperature weak-driving adiabatic 
regime, when n < T, we showed^ that T ; e / / = T lP . 
Then, the conclusion of our previous investigations is 
that for the weak driving adiabatic regime the effective 
temperature defined from a single-particle FDR coincides 
with that determined by a thermometer. 



A non-equilibrium fluctuation-dissipation 
relation 



We analyze the role of effective temperatures (T e ^*) 
from a FDR in the framework of two-particle correlation 
functions. As a priori they are not necessarily the same 
as the effective temperatures defined above we use an as- 
terisk to refer to them. We will focus on current-current 
correlation functions since they are more easily accessible 
from an experimental point of view. We are particularly 
interested in a local relation, that is both currents eval- 
uated at the same point. 

As in the case of the single-particle FDR we focus on 
the dc components of the correlation functions to define 
the effective temperature. This corresponds to assuming 
that an equilibrium-like FDR holds for the k = Floquet 
component with f3 e ff* playing the role of the inverse of 
temperature, 



0ta w 



C*(0u,)-il — Cl(0V) tanh 

(19) 

An equivalent expression for the FDR given in Eq. 
(|19p is obtained by considering the imaginary part, which 
leads to 



tC r £,(0,w) = cath 



Q e f f* 



KM, 



(20) 



where Tp* a (io) — — 21m [C^ a (0, a;)] . (Notice that the real 
part is simply derived by means of Kramers-Kronig re- 
lations.) As in the case of the single-particle FDR in 
Eq. (Tl8| . Eq. (|20)) defines an effective temperature that 
might depend on uj, so the limit u> — > is taken. 

It is important to notice the similarity of this expres- 
sion with the one shown in Eq. (|18l) for single-particle 
Green's functions (fermionic operators). In this case the 
hyperbolic tangent is replaced by an hyperbolic cotan- 
gent due to the bosonic statistic of current operators. 



B. Current-current correlation and noise 

Although we are more interested in the case of local 
current correlations, let us start by considering the more 
general case of correlation at different points. If we con- 
sider two reservoirs (a and /3) and two times (an absolute 
time t and a relative time r) we can define the correlation 
function of currents as 

P a p(t,t-r) = -(AJ a {t)AJ (t-r)+AJ p (t-T)AJ a (t)}, 



(21) 



where AJ a (t) = J a (t) - (J a (t)}. 
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With the definition of the contour-ordered current- 
current correlation function given in Eq. ((5]), the cor- 
relation function of currents given in Eq. (f2Tj) can be 
expressed as 



Pefi&t-T) = \ {C> p {t,t-T)+C<^t-T)) 



(22) 



If instead of a symmetrized current-current correlation 
we are interested in a nonsymmetrized one, 

P$(t,t - t) = (AJ a (t)AMt - r)>, (23) 

the correlation becomes 



P™ {t,t-T)=iC< a {t-T,t) 



In this work we will give results for the symmetrized cor- 
relation only but it is straightforward to obtain the re- 
sults for the nonsymmetrized one. 

Since experimentally the noise spectrum is averaged 
over the absolute time t, the relevant quantity here is 



V a p{uj) = 2 / dT(P a/3 (t,t-T)) t e l 



(25) 



where (. . .)t denotes the time average. From the defini- 
tion of the Floquet-Fourier components given in Eq. (|8]) 
it is easy to see that 



V aP (u) = iC2p(0,u>) 



(26) 



Hence, the only relevant Floquet-Fourier component is 
the one with k = 0. 

As we are considering non-interacting electrons the 
contour-ordered propagator given in Eq. ^ can be ex- 
actly evaluated in terms of single-particle propagators 
((3]). Using Wick's theorem (see the Appendix for the 
details), this contour-ordered function can be written as 



iC a p(t,t') = -W ca W c i3 2_j {Glf3,ka(t' ' ,t)G 'la,kp(t,t') 
ka,k/3 



— Gk/3,ka(t', t)Gi a jp(t, t') 

— Gipjait', t)Gka,kp (t, t) 
+ Gkf3,la(t', t)Gka,lf3(t, t')} 



(27) 



In the Appendix we show the detailed calculation lead- 
ing from this expression to the Floquet-Fourier compo- 
nents (7^(0,0;). Here we only reproduce the results for 
two cases of particular interest. 

The first case is the zero- frequency limit of (7^(0, a;), 
which reads 

P a/3 = l -C%(h = 0,u = 0) = 5 afj P a + P a ^ p , (28) 
where 

/OO J I 00 
=-r a (u') £ ]Tr 7 K)/ Q >',a4) 

■°° fc=-oo 7 



x|Gf Q , 7 (-fc,^.)| 2 , 

L J -°° Z7T k=-oo ^ 

xRe[G^ la (k,u')G^(-k,u}' k ) 

OO 

k— — oo 7 

xGf ath (-k',u' k ,)G^ h (k-k',u' k ,r 

+G> M (k,u')G< M (k^r} 

+ |same with a o (29) 



(24) being 



j afS {u,J) = / a H(i-/„M)+/^M(i-/ a H). (30) 

It is important to notice that this is a general result 
for multiterminal quantum driven systems and the sum 
over 7 extends over all reservoirs connected to the central 
system. For this work we chose a two terminal system, 
but this result is completely general as no assumption 
concerning the reservoirs was made in the calculation. 

At this point it is interesting to compare with previous 
results obtained within the scattering matrix formalism 
(see Ref. [29]). In order to do so, we need to assume 
that all reservoirs are at equal temperature and chemical 
potential (unbiased pump). We split the zero frequency 
noise into two contributions 



__ p (th) j_ p(sh) 

a/3 ' 



P afS EE P^> + P 



(31) 



where 



(sh) 
a/3 



/oo j I 00 
_/( w ')(i-/( w '))r(o/) r (^) 

~°° k— — oo 

Sap £ (\Gt h (k, + \G^ h (-k,cj' k ) 

7 

-|G£, I/3 (fc,a/)| 2 -|G ; g ila (fe,a;')| a } ) 

/OO 1 I 00 ( 

x|Gf Q , 7 (A : , w ')| 2 -/(^) 2 (|Gf C[ ^(fc, W ')| 2 



+ \G%, la (k,J)? +2/(u/)/K)Re Gf Wa (fc,u/) 



- 2Rc 



(/(a/)Gw a (*V) 



-/K)G^(-fe ! ^)*)G< )^a (fc,a;0 , 



-|G<«, Q (fc, w ')| 2 



(32) 



The first term P^p is the Nyquist-Johnson noise while 
Paff 1 IS t nc sn °t n °ise. Using the relation between the 
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Floquet S-matrix and Green's functions 32 

SF,ap(u m , ui n ) — S a p5 n<m — iy/T(uj m )T(uj n ) 

xG£ itf (m-n,w n ), (33) 

it is easy to show that the result given in Eq. (|32p coin- 
cides with the one obtained using the Floquet S-matrix 
formalism in Ref. [29|. 

The other case of interest is the one in which a = 
j5 = P, i.e. we concentrate in current fluctuations of the 
probe. Using the fact that the probe is noninvasive, we 
only keep terms to the lowest order in the coupling w c p 
between the system and the thermometer, 

% E E r >'){/>') 

"°° k=-oo~j=L,R 

x[2-/p(a4+ w )-/ P (u/ fc - w )] 
+ [/pK+ W ) + /pK-c)] 

x(l-/>0)}|G&, /7 (*,a/)| a . (34) 
On the other hand we need Tp* p (uj), which is 

r P ( U ) = r P r^-± 2 r >') 

J ~ oa fe=-oo7=L,fl 

x[/p(wfc-w) - fp(uj' k + U})] 
x\Gf P ^(k,co')\ 2 . (35) 

The functions entering Eqs. (|34p and (|35p are the ones 
involved in the definition of effective temperature given 
in Eq. (gDJ). 

V. RESULTS 

In this section we present results for a central de- 
vice consisting of non-interacting electrons in a one- 
dimensional lattice: 

H = ~wJ2(4^+H.c), (36) 
l,l> 

where w denotes a hopping matrix element between 
neighboring positions 1,1' on the lattice. The driving 
term is chosen as 

2 

Hv(t) = '£eV j (t)cl j c lj , (37) 

with Vj(t) = Eb + Vocos(r2o^ + Sj), Ij being the posi- 
tions where two oscillating fields with frequency Qq and 
phase-lag S are applied. This defines a simple model for 
a quantum pump where two ac gate voltages are applied 
at the walls of a quantum dot i 20 ' 32 i 37 



A. Equivalence between effective and local 
temperature at weak driving 

As in Refs. [l^ and [23| we are interested in the weak 
driving regime, which corresponds to a situation where 
the ac voltage amplitudes are lower than the kinetic en- 
ergy of the electrons in the structure and the driving fre- 
quency is much smaller than the inverse of the dwell time 
of these electrons. We have shown that in this regime the 
local temperature defined from Eq. (|12p. with the chem- 
ical potential of the probe fixed, is identical to the local 
temperature defined from Eq. (|16l) . where the chemi- 
cal potential of the probe has to be determined in order 
to satisfy both equations, and it is also identical to an 
effective temperature defined from a local fluctuation- 
dissipation relation of single-particle Green's functions 
(see Eq. 0). 

We now turn our attention to the effective tempera- 
ture T e ff* defined in Eq. (|20[) . involving current-current 
correlation functions. The correlation functions given in 
Eqs. and (|35l) depend on the temperature Tp and 

the chemical potential pip oi the probe via the Fermi 
function fp. Thus, the effective temperature T e ^*, so 
calculated, also depends on Tp and /ip. There are many 
possible reasonable choices for the latter quantities. In 
this subsection we will concentrate in only one choice and 
leave for the next subsection the analysis of other possi- 
bilities. We choose fip equal to the chemical potential /i 
of the reservoirs and Tp equal to the local temperature 
Tip, i.e. the one for which the heat flow between the 
system and the probe vanishes. 

In Fig. [5]we show a typical plot for iC PP (0,u), Tp* P {uj) 
and their ratio as a function of uj. According to the def- 
inition of effective temperature given in Eq. (f20l) . the 
derivative of this ratio at oj = corresponds to (3 e ^* /2. 
This derivative is calculated numerically. In the same fig- 
ure we plot tanh [(3 e ^*oj/2\ and we see that the quotient 
Tp P (uj)/iC PP {Q,Lu) is well fitted by a FDR-type relation 
for a reasonably large frequency interval. 

In Fig. [3] we show the behavior of the effective tem- 
perature T e M* and the local temperature Tip, calculated 
for the site connected to the left reservoir, as a function 
of the driving frequency fio for two different tempera- 
tures of the reservoirs. This analysis can be done for 
any site of the central system but we chose this particu- 
lar site because its local temperature determine the heat 
current that flows into the left reservoir^ Results for 
any other site of the central system are similar. In Fig. 
[3] the upper panel corresponds to T = 0.016, while the 
lower corresponds to T = 0.005. We see that both ways 
of defining the temperature coincide at low frequencies. 
This supports the idea that, for a given temperature T 
of the reservoirs, T e ^* is a bona fide temperature within 
the low driving regime. As we can see from Fig. [31 the 
higher the temperature T of the reservoirs, the broader 
the region of low driving frequency Qq in which the two 
definitions of the temperature agree. 

As we mentioned earlier, the definition given in Eq. 
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FIG. 2: (Color online) Current-current correlation functions 
Tp* p {uj) (dotted red), iCp P (0, u) (dashed blue), their quotient 
(green diamonds), and tanh[/3 e -"*a;/2] (solid black) as a func- 
tion of lu. The reservoirs have chemical potential fi — 0.2 and 
temperature T = 0.025. The driving frequency is fio = 0.01, 
the amplitude is Vo = 0.05 and Eb = 0.2. 



t i i | i i i | i i i | i i i | i i i | i i r 




FIG. 3: (Color online) Local temperature Tip (dashed black) 
and effective temperature Tip (solid red) for the site IP = 
IL (i.e. the site connected to the left reservoir) as a function of 
driving frequency flo- The reservoirs have chemical potential 
fj, = 0.2. The upper panel corresponds to T — 0.016, while 
the lower panel corresponds to T — 0.005. 

([20)) can be used to calculate the effective temperature 
in each site of the central system. In Fig. 2] we show the 
comparison between the local temperature Tip and the 
effective temperature T e ^* all along the sample. The 
values of Tip and T e ^* are plotted for each point of the 
linear chain, for Tl = Tp = T = 0.02, = pip — 
fi = 0.2 and a particular low value of Oo = 0.001. We 
can see that there is a good agreement between the two 
temperatures along the whole structure and an almost 
perfect agreement within the "Left" and "Right" regions 
(defined in Fig. QJ , which are the ones from where we can 
determine the heat flow between the system and each one 
of the reservoirs (see Refs. I23lf33f ). It is also important 
to notice the existence of 2kp Friedel-like oscillations, Uf 
being the Fermi vector of the electrons leaving the reser- 




FIG. 4: (Color online) Local temperature Tip (black dia- 
monds) and effective temperature T^p 1 * (red circles) along a 
one-dimensional model of N — 30 sites with two ac fields op- 
erating with a phase lag of S = n/2 at the positions indicated 
by dotted lines. The system is in contact with reservoirs with 
chemical potentials fi = 0.2 and temperature T = 0.02. The 
driving frequency is Qq = 0.001, the amplitude is Vo = 0.05 
and E B = 0.2. 



voirs. These oscillations are an indication of quantum 
interference. They were previously reported for exactly 
the same setup we study in this wor h 19 i 23 and also pre- 
dicted in other mesoscopic systems under a stationary 
driving. 38 



B. Different choices of Tp and fip 

The effective temperature T e ff* depends on the values 
of Tp and fip (respectively the temperature and chemical 
potential of the probe). The choice analyzed in the pre- 
vious section was [ip equal to the chemical potential /x of 
the reservoirs and Tp equal to the local temperature Tip. 
We will call this choice Case I. Another suitable choice 
(Case II) could be to choose /ip = (i, as in the previous 
case but T P such that T eff * = T P . A third choice (Case 
III) could be to choose T P such that T eff * = T P but at 
the same time /ip = flip (the local voltage) in order to 
have a vanishing charge current between the system and 
the probe at that temperature. In this work we will only 
deal with these three possibilities. 

If Fig. [5] we show the three different effective tem- 
peratures corresponding to the above mentioned cases 
together with the local temperature as a function of the 
driving frequency Qo for a given temperature of the reser- 
voirs. As we can see, all three cases give a good estimate 
of the local temperature in the regime of interest (i.e. 
low driving frequencies). This behavior supports the ro- 
bustness of the definition of the local temperature from 
a FDR. 
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FIG. 5: (Color online) Local temperature Tip (dashed black) 
and effective temperature T ; p for Case I (solid red), Case 
II (dotted blue) and Case III (dashed and dotted orange) for 
the site IP — IL as a function of driving frequency Qo. The 
reservoirs have chemical potential [i — 0.2 and temperature 

r = o.oi. 



VI. SUMMARY AND CONCLUSIONS 

In this work we have calculated the current-current 
correlation functions for quantum driven systems and 
found an explicit expression for the zero-frequency noise 
within the Schwinger-Keldysh Green's functions formal- 
ism. In the particular case of multiterminal unbiased 
quantum driven systems our result is in agreement with 
previous results obtained within the scattering matrix 
approach.— For non-interacting systems both descrip- 
tions agree, while the Green's functions has the advan- 
tage of providing a systematic framework for the study 
of interacting systems. 

We have also defined an effective temperature from a 
local fluctuation-dissipation relation for current-current 
correlation functions and showed that for low frequen- 
cies it coincides with the local temperature defined with 
a thermometer and from a FDR at the level of single- 
particle propagators. This result opens the possibility of 
using current-current correlation in real experiments, in 
order to define the local temperature of a driven sample. 
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Appendix: Analytical expressions for 
current-current correlation functions 



Using Eq. (0 , the connected contour-ordered current- 
current correlation (see Eq. ([5])) can be written in terms 



of electron operators as 



—W ca W c p £ (^D ka ,1a, k 13,1/3 (*)*') 
ka,k(3 

— Dka,la,l/3,kl3(t, t) 



+Di ai ka,l/3,kp(t, t')j > 



(A.I) 



where 



D itj , k ,i(t,t') = (Tciclit^clit'Mt')]) 

Using Wick's theorem and the definition of the 
contour-ordered Green's function (Eq. ([3])) we can 
rewrite Eq. (|A.2|) as 

D itj , k ,i(t,t') = G hl {t',t)G^ k {t,t'). (A.3) 

The substitution of Eq. (|A.3j) into Eq. (|A.L|) gives the 



result shown in Eq. (|27|) . 

The next step is to calculate the lesser and greater 
Green's functions. In order to have the lesser (greater) 
Green's function we need t E G\ and t 1 £ C2 (t € C2 and 
t' e Ci). If t £ Ci and t' e C 2 , then 

D<j M W) = G U*> t ) G H t > i! )> (A.4) 

Using Eq. (IA.4[) we can write the lesser Green's func- 
tion as 



ka 7 k(3 

_ (Gh atk p{t,t')*Gf a ^ (t,t')J 

-(GL w (t,tTG< a ^(t,t')) 

+ (G? aM (t,trG>: a , l0 (t,t'))}. (A.5) 

where the property Gf'j?(t,t') = -Gf,'^ (t' ,t)* was used. 
Notice that the last term is equal to the first interchang- 
ing < with > and conjugating. The same situation arises 
with the second and third terms. The result for the 
greater Green's function is obtained by switching < with 
>. 

We are interested in the Floquet-Fourier components 
of this Green's functions. So, according to Eq. ([8]) we 
need to calculate 

Substituting with the Floquet-Fourier expansions of 
the Green's functions (see Eq. (|8])) into Eq. (|A.5[) and 
rewriting into the form of Eq. (IA.6|) we can obtain the 
k = Floquet-Fourier component, which is 

C< p (0, u) = i [a< p (u) + A> (-uT + B< p {lo) + B> (-u>y 

(A.7) 
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where 

00 POO (J^jJ 1 

AZpip) = E / ^-E^W-^) 

fe/3 

W = E / ^G? a> w(k,u'y 

k=-oo J -°° 

and A' > is obtained from X < by switching < with > (X = 
A, B) . Since the greater Floquet-Fourier component is 
obtained from Eq. (|A.5I) by switching < with >, it can 
be written as: 

C>p(0,uj) = l [A> al M + A< p (~u)* + B> (u) + B< p (-uy 

(A.8) 

where X > is obtained from X < by switching < with > 
(X = A,B,C,V). 

We want to obtain an expression of the previous quan- 
tities in terms of Gf- V (fc, ui) with I, I 1 in the central system. 
In order to do so we need to calculate the Floquet-Fourier 
components of 



ka 



and 



where 



+G< M (t,t 1 )^(t 1 ,t') , (A.9) 



— - — {5 ap -E<{t,t') + [ dhdt 2 

WcaWc/3 L J 

x[Ef(t,t 1 )G^(ti,t 2 )S<(t 2 ,t') 

(A.10) 



An 



E*(f,f) = Ej(f,t)* 



(A.ll) 



with their respective Fourier transforms 



dw' r a (u/) 



Using the Fourier transforms of E< and E fl together 
with the Floquet-Fourier expansion for the Green's func- 
tions (see Eq. jSJ) we obtain 



E" 



G$,*„(*.w) = ^-{Gg )ta (fc,o;)E<( W ) 

(fc, W )E*( W )*},(A.13) 



where 

00 

m— — oo 7 

xGg, l7 (*-m,w^), (A.14) 
with E 7 (cj) = r 7 (oj)A^(aj), and 



•^H = »/y( W )> 



(A.15) 



On the other hand, 



— (<Ja^4oS^(w) 



U! CQ ,W c/ 3 

+G,*, /J (fc,w)E|(w)E»(w fc ) 

+ [Gg !iQ (-fc, Wfc )*s<K) 

+G<^(fc,a;)E«( Wfc )]sf(a;)*}. 

(A.16) 



By substituting Eqs. (|A~T3|) and (|A~T6)l into Eq. (|A~8|) 
and then into Eq. (|A.7[) we obtain an expression for 
C^g(0, ui). The expression for the greater Green's func- 
tion is obtained by switching < with >. 

We are interested in the case where a — j3 = P, i.e. we 
are interested in fluctuations in the current flowing to the 
probe. As we are dealing with a non-invasive probe, we 
only keep terms up to the lowest order in the coupling be- 
tween the system and the probe. Since T a ,H a oc \w ca \ 2 
it is easy to see that 



A< p (oj) = 0(Kp| 4 ), 

/GO J J 
— G> >iP (0,a/)*A<(u/ + u>)T P {u,' + u>) 

+0(Kp| 4 ). (A.17) 

Thus, Bp P (uj) and Bp P (uj) are the only terms with con- 
tributions of order |w c p| 2 . Using Eqs. (IA.14[) and (|A.15|) 
we can write, up to this order, 



^~ E E \G? P , h {-n,J n )? 

-OO Z ~ ,._ r r> 



Cg P (0,u)) = -i 



k— — 00 7— L,R 



(A.12) 



x {/ 7 p«, w' + w)r P (w' + w) 

+/iP( w l. w '- w ) r F( w '- w )} r 7( w l)- 

(A.18) 



10 



To obtain an expression for ip* P (oj) we can use the fol- 
lowing identity, easily deduced from Eq. (jlOp . 



lp* P (cj) = -2Im[C£p(0,w)] = i (C> P (0,u) ~ C< P (0,w)) , 

(A.19) 

which leads to 



^p(w) = 



E E r 7 K)|Gf F>h (-n,^)| 

— 00 7 — L , i? 



x { - /p(w' + u)] r F (u/ + w ) 

- / 7 K)-/p(w'- w )]r P ( w '- w )}. (A.20) 

If we additionally consider a constant density of states 
in the probe to obtain a result independent of any par- 
ticular probe, we arrive to the expressions given in Eqs. 
(ES and (ESI). 
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